Extraordinary wetting phase diagram for mixtures of Bose-Einstein condensates 
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The possibility of wetting phase transitions in Bose-Einstein condensed gases is predicted on the 
basis of Gross-Pitaevskii theory. The surface of a binary mixture of Bose-Einstein condensates 
can undergo a first-order wetting phase transition upon varying the interparticle interactions, using, 
e.g., Feshbach resonances. Interesting ultralow-temperature effects shape the wetting phase diagram. 
The prewetting transition is, contrary to general expectations, not of first order but critical, and the 
prewetting line does not meet the bulk phase coexistence line tangentially. Experimental verification 
of these extraordinary results is called for, especially now that it has become possible, using optical 
methods, to realize a planar "hard wall" boundary for the condensates. 

PACS numbers: 03.75.Hh, 68.03.Cd, 68.08.Bc 
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In this Letter we pose and answer theoretically the fol- 
lowing fundamental questions. Is a wetting phase tran- 
sition possible in ultra-cold dilute gases which undergo 
Bose-Einstein condensation (BEC)? How many species 
of atoms are needed to compose two coexisting phases 
with an interface? What is the nature of the surface or 
"wall" at which these phases are "adsorbed" ? If a wetting 
transition occurs, is its character first-order or critical? 
What is the nature of possible prewetting phenomena 
away from bulk two-phase coexistence? We focus on the 
essential physics in the application of wetting theory to 
BEC and give details of the calculations elsewhere^. 

The simplest system, a one-component gas, is insuffi- 
cient for studying wetting transitions, because the con- 
densate fraction and the normal fraction are fully mixed 
in position space. Spatial segregation is only possible in 
an external potential (e.g., gravity) 2 . Without this po- 
tential an interface between normal fraction and conden- 
sate does not exist, and therefore the essential interfacial 
tension is missing. 

The next simplest BEC system is a two-component 
gas. Binary mixtures of trapped BE condensates 
of alkali-metal atoms have received much attention, 
experimentally 3 ^ and theoretical!^^ 7 ", since the seminal 
predictions concerning their phase behavior by Ho and 
Shenoy- . Different degrees of spatial segregation of two 
pure-component condensates at two-phase coexistence, 
and the possibility of a symmetric-asymmetric configu- 
rational transition in a trap, have been elucidated by Ao 
and Chui 6 and Svidzinsky and Chui 7 . 

We exploit these findings in the context of wetting 
phenomena&2ii2iii, and point out the existence of a sur- 
face phase transition from partial to complete wetting, 
upon varying the intra- or inter-species interactions in 
mixtures of BE condensates. The surface consists of a 
planar "hard wall" , at which the condensate wave func- 
tions vanish. 

The basic physics of wetting is best indicated by invok- 
ing the familiar energy balance known as Young's law, 



where "f wi is the surface (free) energy of a phase of pure 
component i, 7 12 is the interfacial tension between pure 
phases 1 and 2, and is the contact angle with which the 
1-2 interface meets the surface. One of the components, 
say 2, is preferentially adsorbed at the surface, so that 
Iwi > 7w2- The surface phase transition from partial 
wetting (0 > 0; Figure la) to complete wetting (9 = 0; 
Figure lb) is then the dramatic phenomenon in which, 
at the surface, pure phase 1 is displaced by a macroscop- 
ically thick wetting layer of pure phase 2. Thus, the Wl 
surface is replaced by a W2 surface in parallel with a 
1-2 interface. The physical implications of the singular- 
ity structure around wetting transitions have been the 
subject of impressive experimental^, theoretical^ and 
simulational 13 research. 
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FIG. 1: a) Partial wetting. The interface between the two 
condensates makes a finite contact angle with the wall, b) 
Complete wetting. A macroscopic layer of condensate 2 in- 
trudes between the wall and condensate 1. 



In view of the ultra-low (typically nanoKelvin) tem- 
perature needed for BEC in dilute gases, the quantum 
mechanical many-body theory at T = is appropri- 
ate. For weakly interacting gases the Bogoliubov mean- 
field theory of BEG 2 *^^ is justified. For our inhomo- 
geneous binary mixture one obtains two coupled Gross- 
Pitaevskii (CP) equations for the condensate order pa- 
rameters V'l = Fie lXl and ■02 = -F^e 1 * 2 . We recall that ip 
is the ground state expectation value of the Boson field 
operator. It can be chosen real here, since stationary 
(non-flowing) condensates are assumed, so \i — X2 = 0. 
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We adopt the standard geometry for wetting and con- 
fine the atoms to the half-space z > 0. The relevant 
surface is the z = plane. It is then natural to employ 
the grand canonical ensemble, with chemical potentials 
pi and p-J^. The GP equations are the Euler-Lagrange 
equations of the grand potential functional 
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The sum is over the two species, with atomic masses 
m,; . The repulsive interactions Gij > are related to the 
s-wave scattering lengths through Ga — AttH au/irii 
for like particles and Gyi = Airfi 2 ^^(l / 'mj + l/m-2) for 
unlike particles. For alkali-metal atoms, a sa 10 2 A. 

The confining potential is taken to be a hard wall at 
z = 0, so that 



^iO,J/,0) = ip 2 ( x ,y,0) = 



(3) 



The closest experimentalists have come to make a hard 
wall is to use an evanescent wave, blue-detuned, extend- 
ing from within a planar prismiS. The turn-on distance 
of this optical wall is typically only about X/2tt « 80nm. 
To contain the atoms, a conventional quadratic-confining 
magnetic trap can be added, so that a "quasi-square" 
exponential potential results for z < 0, and the usual 
ax 2 + by 2 + cz 2 for z > 0. In order for volume forces 
due to a non-uniform external potential to be negligible 
compared to surface forces governing wetting, the har- 
monic potential must be sufficiently flat-bottomed near 
r = 0. For a characteristic length of the order of microns 
or longer we can ignore the harmonic potential in the 
calculations and assume translational invariance of ip in 
x and y directions. 

The known condition for bulk phase separation of the 
binary mixture is that the unlike particles repel each 
other more strongly than the like ones on average, K = 
G12/VG11G22 > 1- Otherwise, a single mixed phase re- 
sults. Bulk phase coexistence requires equal pressure, 
Pi = P2, in the two phases, with Pi = pi/2Gu. A bulk 
triple point with coexisting mixed and pure phases 1 and 
2 is found for K = 1. This triple point bears a remi- 
niscence to a critical point in that the interfacial tension 
between pure phases 1 and 2 vanishes. However, the 
phases themselves remain distinct. 

The bulk condensate number densities in the pure 
phases are poi = I^M 00 )! 2 = HlGa. A "healing length" 
£i = hj y/2mipoiGu characterizes the recovery distance 
of the order parameter from a disturbance, e.g., a sur- 
face where p = 0. A typical value is £ ~ 10 3 A. The 
two surface grand potentials j wi are proportional to 
the product of the pressure and this healing length^, 
Iwi = (4\/2/3)P£i, which defines a surface "thickness". 



For strongly repulsive interactions and/or high densities, 
the healing length of a condensate is short, and its surface 
tension is low, so that it prefers to be near the surface. 
This suggests to define a "surface field" proportional to 
the difference £1 — £2 • 

The second important length is the "penetration 
depth" A, characterizing the distance over which one con- 
densate decays or "penetrates" into the other at the 1-2 
interface^, Aj = £i/y/K — 1. A typical value is A ss 1pm. 
This length diverges when approaching the triple point. 
Ao and Chui discuss the following two limits for which 
we determine the wetting behavior. 

A) Strong segregation limit: Partial Wetting 
In the limit A -C £ or K — *■ 00 the condensates show 
almost no mutual penetration. The interface then con- 
sists of two surfaces, one where condensate 1 decays to 
zero over the healing length £1, and an adjacent one 
where condensate 2 vanishes over £2 ■ The interface thick- 
ness is thus £1 + £2 and its tension is approximately 
7i2 ~ 7wi + 7w2- There can be no complete wetting, 
since the energy cost of an interface is too high to be 
able to satisfy Young's law with 8 — 0. 

B) Weak segregation limit: Complete Wetting 
In the experimentally more relevant limit A > ^ or 
K — ► 1 the mutual penetration of the condensates leads 
to mixing on the scale of Ai + A2. In spite of its huge 
thickness, the interface has a low tension, since the en- 
ergy per unit volume scales asi& (£/A) 2 P. As a result, 
7 12 oc + &)VK — 1- The interfacial tension van- 
ishes at the triple point with a square-root singularity. 
For Young's contact angle we obtain 



cos 9 = const 



(4) 



This implies that a wetting transition is unavoidable 
upon approach of the triple point. Indeed, the wetting 
phase boundary implied by J3J is given by 



(£i/6-i)«(^-i) 1/2 



(5) 



A full numerical computation^ corroborates the phys- 
ical insights gained so far, and leads to the wetting 
phase diagram at bulk two-phase coexistence shown in 
Fig. 2. The wetting transition is of first order. On the 
phase boundary the grand potential of partial wetting 
crosses that of complete wetting, and both states have 
metastable continuations. Interestingly, at the wetting 
transition the grand potential is degenerate: all wetting 
layers have the same energy, regardless of the layer thick- 
ness. This is akin to what happens at the wetting transi- 
tion in an Ising model at T — 0, where the interface con- 
sists of a plane of broken bonds at an arbitrary distance 
from the surface. The usual entropic repulsion between 
the surface and the interface due to interfacial capillary 
wave fluctuations is absent at T = 0. 

Yet, in contrast with Ising spins on a lattice, the con- 
tinuum field theory we are dealing with allows one to 
"inflate" a wetting layer by varying its shape and thick- 
ness continuously between infinitesimal and macroscopic, 
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FIG. 2: Wetting phase diagram in the plane of surface field 
(£i — ^2)/^2 and relative interaction strength K. At bulk 
two-phase coexistence (K > 1) a first-order phase boundary 
separates partial wetting from complete wetting, for £i/f2 > 
1. For £1/^2 < 1 the roles of the condendates are interchanged 
and one may use "drying" in place of "wetting". The phase 
boundary is parabolic near the triple point (1, 1). 

at constant grand potential. A remarkable consequence 
of this is the coincidence, at bulk two-phase coexistence, 
of nucleation and wetting transitions^. For illustration, 
Figure 3 shows two order parameter profiles for parame- 
ters K and ^1/^2 on the wetting line. 




FIG. 3: Inflation of a wetting film of condensate 2 at the 
surface of condensate 1. Normalized order parameters, fi = 
Fi/y/poi, are shown for two surface states with the same grand 
potential, at first-order wetting (K — 1.5, £1/^2 = 2). 

A crucial physical question is how the surface field 
and/or the relative interaction strength can be tuned 
experimentally. Importantly, these variables can be ex- 
pressed in terms of the atomic masses and scattering 
lengths alone, and are therefore fully microscopic. Fesh- 
bach resonances allow manipulation of the scattering 
length(s) over more than an order of magnitude^ and 
are a suitable tool for exploring the phase diagram. 

Away from bulk two-phase coexistence, a first-order 
wetting transition normally possesses an extension into 
the bulk one-phase region where the wetting phase is 
metastable. If this extension exists, general arguments 
show that it must also be of first-order, at least from 
the wetting point up till a prewetting (multi-)critical 
point in the phase diagram. Furthermore, the prewet- 
ting line must meet the bulk two-phase coexistence line 
tangentially 20 . To our knowledge, in all systems that 
have been observed or predicted to undergo a first- 
order wetting transition, ranging from classical adsorbed 



fluidsii or superfluidsi^, over magnetic systems^, to 
type-I superconductors^, the prewetting transition is 
also of first order. Yet, mixtures of BE condensates pro- 
vide us with an exception, as we now show. 

Figure 4 shows the prewetting phase diagram in the 
space of "bulk field" and relative interaction strength K. 
The bulk field which can drive the mixture away from 
bulk two-phase coexistence is proportional to the pres- 
sure difference P\ — P 2 . The combination \fP\JP~2 — 1 
serves the same purpose. Note that pure phase 1 (2) is 
stable in bulk for P\ > (<)^2 and that the mixed phase is 
stable for K < mm{ y /P 1 /P 2 , y/P 2 /Pi}. The bulk triple 
point is at (1, 1). For a given value of the surface field, 
i.e., for fixed the wetting point W and the prewet- 

ting line are shown. When, along the path indicated by 
arrows, the prewetting transition is crossed, an infinites- 
imal film of condensate 2 is nucleated at the surface of 
condensate 1. Moving further to the left, the thickness 
of this film grows and, upon approach of two-phase co- 
existence (at Pi = P2), it diverges logarithmically. This 
slow divergence is expected for the approach to complete 
wetting in systems with short-range interactions. The ex- 
traordinary physics emanating from this phase diagram 
concerns i) the second-order or "critical" nature of the en- 
tire prewetting line, and ii) the fact that the prewetting 
line meets bulk two-phase coexistence at a finite angle 
(clearly seen from the mathematical continuation along 
the dashed line). 



6 \ 




FIG. 4: The phase diagram for prewetting at £1/^2 = 4, 
in the space of relative interaction strength K and bulk field 
\fp7JP2-l. The bulk phases " pure 1" , "pure 2" and "mixed" 
coexist at the triple point T. Bulk two-phase coexistence of 
the pure condensates is along the line Pi = P2, on which 
the first-order wetting transition W is indicated. The second- 
order prewetting phase boundary intersects the bulk transi- 
tion at W . The line with arrows indicates a possible approach 
to complete wetting. 

If there is nothing wrong with the general thermo- 
dynamic arguments concerning prewetting, how is this 
anomaly possible? There are two premises in the general 
reasoning, which are not fulfilled here. Firstly, at this 
first-order wetting transition there is no abrupt jump of 
the film thickness between a microscopic and a macro- 
scopic value, but all intermediate values are equally sta- 
ble due to the grand potential degeneracy. Since there is 
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no interface potential barrier, for consistency the prewet- 
ting transition must be continuous instead of first-order, 
and we find that it is. Secondly, since this prewetting 
transition is a second-order nucleation transition, unlike 
for ordinary prewetting, the difference in adsorption on 
either side of the transition does not diverge upon ap- 
proach of the wetting transition at W, and the condition 
for a tangential meeting is not met. Consistently, we find 
an intersection at a finite angle. 

Within the wetting phenomenology we reinterpret the 
interesting symmetric-asymmetric (SA) transitions pre- 
dicted for an axially symmetric square-well trap&, and 
for a conventional quadratic magnetic trapi In the sym- 
metric (S) state the condensates form concentric clouds. 
This corresponds to complete wetting (Fig. lb), with, e.g., 
a W2 surface and an embedded 1-2 interface in a "cherry" 
configuration. This state may exchange stability with 
one (A) which breaks the trap symmetry, with a planar 
1-2 interface cutting across the system. This is partial 
wetting (Fig. la). 

For a square- well trap, a surface energy balance deter- 
mines the SA transition 6 . However, it is not governed 
by Young's law, but by finite-size effects, cylindrical ge- 
ometry and particle number conservation. Consequently, 
the contact angle at the SA transition does not tend to 
zero continuously as in the bona fide wetting transition. 
Instead, it jumps from a finite value 0$a to zero. For a 



quadratic trap, surface and volume contributions are en- 
tangled in the energy balance^. In the S state the com- 
ponent with the largest self-repulsion is on the outside. 
For sufficiently small trap radius, an A state is favored, 
to relieve the high capillary pressure across the curved 
1-2 interface in the S state. The wetting interpretation 
is complicated here due to the spatially varying external 
potential and the absence of an articulated surface. 

In conclusion, in view of the extraordinary ultralow- 
temperature effects in the predicted wetting phase di- 
agram for binary mixtures of Bose-Einstcin condensed 
gases, we advocate an experimental study of wetting and 
prewetting in a trap which is suitable for observing true 
wetting singularities. The confinement should consist of 
a half-space with one planar "hard wall", i.e., a steep 
repulsion with a turn-on length smaller than the heal- 
ing lengths and penetration depths of the condensates. 
A mild conventional trap may be used to keep the gas 
near the wall. The wetting transitions can be induced by 
tuning the interatomic scattering lengths using Feshbach 
resonances. 
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